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Abstract 

Electrohydrodynamic (EHD) generation, a commonly used method in 
BioMEMS, plays a significant role in the pulsed-release drug delivery 
system for a decade. In this paper, an EHD based drug delivery system is 
well designed, which can be used to generate a single drug droplet as small 
as 2.83tjX in 8.5ms with a total device of 2 X2x3mm 3 , and an external 
supplied voltage of 1500V. Theoretically we derive the expressions for 
the size and the formation time of a droplet generated by EHD method, 
while taking into account the drug supply rate, properties of liquid, gap 
between electrodes, nozzle size, and charged droplet neutralization. This 
work proves a repeatable, stable and controllable droplet generation and 
delivery system based on EHD method. 
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a, [m] Radius of nozzle 

d, [m] Gap between tip of nozzle and ground electrode 
E, [V/m] Electric field strength 

F e , [N] Electric force 

F s , [N] Surface force 

h, [m] Height of liquid droplet cap 

m,[g/s] Drug supply rate 

r, [m] Radius of droplet pending at tip of nozzle 

r e f j , [m] Effective radius of droplet 

U, [V] Applied external voltage 

u, [m/s] Velocity of drug fluid in nozzle 

7, [N/m] Surface tension 

e, [F/m] Drug fluid permittivity 
9 C , [rad] Contact angle 

/i, [Pa • s] Drug fluid viscosity 

£, [V] Surface electric potential 

p e ,[C/m 3 ] Spatial charge density 
Pd, [kg/m s ] Drug density 

1 Introduction 

Microfluidics and related nano/micro-technology have been active research 
fields for several years. Numerous microflow devices such as micro channels 
[HE], flow sensors PQEBUT], microvalves [3], micropumps ffi [5].[9]-[T2]. 
electrowetting [H [6], electro-hydrodynamics (EHD) [13]-[15j,[21J-[24] and 
magnetic-hydrodynamic [7] have been reported. A droplet dispenser using 
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the microfluidic chip format [8] has recently found increasing popularity 
due to its potential for massively parallel array ability and scalability of its 
size |18| . Especially, the latter characteristic reduces the dead volume of 
the device, which is essential when dealing with an extremely small sample 
volume. 

The pulsed and discontinuous drug delivery systems are more advanta- 
geous in the academic and clinic applications, which can be used in ther- 
apeutics as concentrations of solutions vary with time, or matches bodys 
release of peptides and hormones [TJ 119] . The device would deliver drugs 
according to various stimuli such as chemical pH, electric fields, and tem- 
perature. By judging the degree of external signal, and the device would 
release appropriate amounts of drug [201 121] . 

In this paper, the EHD based drug delivery system is applied mainly for 
three reasons. Firstly, it doesnt require complex fabrication because there is 
no moving part in the whole device. Secondly, it can produce mono-disperse 
droplets in a wide range of sizes using either electrospray or single-droplet 
generation. Thirdly, instability, poor directionality, and diversity in droplet 
size can be reduced by coating Teflon on surface of PDMS chip and reducing 
surface tension of liquid in fabrication process |21[ 123] . Our work shows that 
a stable, repeatable and controllable generation of volume at a constant rate 
has been achieved by the EHD method. Many types of solutions can be 
delivered through the device such as DI water, acetic acid, methanol, and 
so on. 



2 Device Concept and Operating Principle 

A typical EHD based micro drug delivery device is shown in Fig. [TJ An 
electric field would be formed after when we apply voltage on two electrodes. 
Drug solutions would be under the electric pressure as soon as they arrive 
at the tip of the red nozzle. Under the pressure, the liquid is elongated and 
breaks up into single micro-droplets from the tip. Taylor cone, which is very 
common in hydrodynamic spray processes, is suppressed by controlling the 
surface wetting property of the PDMS device and the surface tension of the 
sample liquids. 
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Figure 1. Schematic of EHD drug delivery device 




Figure 2. Design and dimensions of EHD based device, (a) 
PDMS core (blue), external box (grey) and applied voltage 
(black); (b) dimensions of PDMS core and metal tube. 



3 Model Design and Fabrication 



3.1 Device Design and Dimensions 



The core of the device made of PDMS consists of a spherical drug reservoir 
and a cylindrical tube extruded from the reservoir, as shown in Fig. [2j A 
metal tube/needle, as an anode, is inserted into the PDMS core. 



3.2 Materials 



Electrode: Silicon-based electrode coated with polyimide over an adhesion 
promoter. 

Insulating material: Diamond-like carbon (DLC) film adheres strongly to 
some plastics providing DLC is biocompatible. 
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Figure 3. Fabrication process of drug reservoir and Teflon 
coating at tip of nozzle 



3.3 Fabrication 



The big chunk of PDMS in the middle is supposed to be made up by two 
half parts. For the reservoir, it is a half sphere, and for the micro-channel, it 
is a half cylinder. Overall, the fabrication will be complicated. So we choose 
to use Grey Scale method on SU-8 of E-beam lithography. Teflon coated 
on the inner and outer surface of tip of the metal nozzle, as shown in Fig. 
[3^d), can improve the stability and directionality of the droplet formation, 
and further reduce the droplet size by reducing surface tension of liquid. 

Fabrication process, as shown in Fig. [3j consists of four steps as follows. 
(1) Put SU-8 E-beam resist on top of silicon wafer. (2) Grey Scale to overex- 
pose the middle, and by manipulating the gradient density to get a perfect 
half cylinder and a half sphere. In Fig. [3^b), 1 and 0.5 [iC/cm 2 , is just 
a trial number. (3) After overexposure, one half core is done. Then cast 
PDMS on top of it. (4) Coat metal nozzle with Teflon. 



4 Charged Droplet Neutralization 



Droplets generated by EHD method would carry a slight amount of charges. 
Since nobody has really taken the droplet neutralization effect into account 
for a matured drug delivery system, this part is hard to accomplish. Here, 
we give out three possible solutions for this neutralization problem. A neu- 
tralization system is shown in Fig. [4j 
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Figure 4. Schematic of a neutralization system 

4.1 Negative Corona Discharge 

When the charged droplet comes out of CEl, by changing HV3, we can 
actually change the negative electrons released from tip N, which are used 
to neutralize the positive charges carried on droplet |25j. 

4.2 Negative Spraying of Volatile Liquid 

By replacing the N part with a capillary of spraying volatile liquid out, which 
will bring negative electrons in the liquid, when the evaporation is finished, 
electrons remain to neutralize the positive charges. 

4.3 Electron Cloud 

Electrons generated at the wall by proton losses are accelerated and decel- 
erated by the beam potential and hit the opposite wall with a net energy 
gain, producing secondary electrons and forming a charged cloud. 

5 Calculation of Droplet Size and Formation Time 
5.1 Droplet Size 

Droplets are formed at the tip of nozzle, that is, the metal tube inserted into 
PDMS channel. Figure [5] below shows the formation process of a droplet. 
Now, assume the geometry of a liquid drop as a spherical cap. It is well 
known that the volume of a liquid droplet formed at the tip is given by 

(1) Vhquid = Vcap = ^h 2 (3r - h 2 ) 
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Figure 5. Formation process of a droplet, (a) (e) 




Figure 6. Parameters used in calculation of a droplet formation 

where, h is height of a droplet in the larger semisphere, h = r + \J r 1 — a 2 , 
and a is radius of nozzle. 

Then, define the spatial charge density, p e , in C/m 3 . The total charge of 
a liquid drop is 

(2) Q = /0 e X Vuguid 

In MEMS device, the gravitational effect on a droplet can be neglected. 
Therefore, there are two groups of forces exerted on a forming droplet, elec- 
tric force (F e ) and surface force (F s ), as shown in Fig. [6j Electric force 
caused by the external electric field, due to Coulombs law, is given by 

(3) F e = E x Q = E Pe V liquU 

where, E = ^, and, U is applied external voltage and d is gap between tip 
of nozzle and ground eletrode. 
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Surface force due to the surface tension between liquid, gas and nozzle, is 
given by 

(4) F s = 27ra7sin(7r - C ) 
where, contact angle, theta c = § + arccos (^). 

A drug droplet is formed at the tip of nozzle, only when the electric force 
caused by EHD effect is larger than the surface force on the meniscus of 
droplet, F e > F s . The droplet is elongated and breaks up from the nozzle 
at force equilibrium. The effective radius of a droplet can be calculated as 
follows. In equilibrium, F e = F s , equating Eqs. (3) and (4), volume of a 
droplet is 

4 3 27ra7sin6» c 

(5) V droplet = ^r £ff = — , 

and effective radius is 

'3a7 sin 9 C \ 



( 6 ) r eff 



2E Pe 



So far, we have calculated the effective radius of a formed droplet in 
equilibrium at static state. Based on our results, the drop size is dependent 
on nozzle size, drug solution properties (surface tension), Teflon coating 
(contact angle), spatial charge density, distance between nozzle and ground 
electrode, and supplied voltage. 



5.2 Droplet Formation Time 



Two different time constants are of our interests. 



tg a p, the required time between two formed droplets, is given by 

,^ ^ pdVdropiet = lia-ypd sin 9 C 

gav m Ep e m 

It significantly depends on the drug supply rate, m. 

i/ormjthe required time to form a single droplet from the tip of nozzle, 

_ Vdroplet _ Vdroplet 



(8) t fc 



f£v(r,0)rdrde J % (r) 2vrrdr 



where, v (r, 6) is fluid velocity profile created at the tip of nozzle. We will 
show how to calculate the velocity profile in the next section. 
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Figure 7. EOF model, (a) drug flow under constant elec- 
tric field; (b) two possible velocity profiles at outlet, parabolic 
(green) and constant (yellow) velocity. 



5.3 Velocity Profile inside Nozzle 



In order to estimate the flow rate inside the electrode induced by an external 
voltage, we use a well established model called Electro-osmatic flow (EOF), 
which is the motion of liquid induced by an applied electric field across a 
charged micro-channel, as shown in Fig. [7j 

Several assumptions are made for the calculation of velocity profile at tip. 
(1) Two-dimensional fully developed flow. (2) Zero pressure gradient. (3) 
Slip velocity at channel wall is u = —E, where, e is drug fluid permittivity, 
\x is drug fluid viscosity, and £ is surface electric potential |26j . 

The governing equations used in our calculation are shown below, 

(9a) V • u = 

(9b) pu-Vu = -Vp + /xV 2 ii 



Write the continuity and momentum equations in Cartesian form, 
du dv 

< 10a > 3i + Wy ~ ° 

, , . f du du\ p ( d 2 u d 2 u 



dx dy J dx \dx 2 dy 2 

, 10 ) , dv ^ dv\ P + f ® 2y _|_ ® 2y 

dx dy J dy \dx 2 dy 2 



Simplify the equations, we have 

(11) v = -»■ ^ = -»■ u = u(y). 

dx 
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Using the assumption of zero pressure gradient, Eq. (10c) vanishes and 
Eq. (10b) can be simplified as: 

d 2 u 
l dy 2 



(12) = M - 



A general solution for Eq. (13) is in a form of 
(13) u = C lV + C 2 

Here, C\ and C 2 can be calculated, due to the boundary conditions, from 

(14a) u{-a) = -—E = -C 1 a + C 2 

ft 

(14b) u{a) = -—E = C 1 a + C 2 

ft 

Therefore, d = and C 2 = -fE. 

The fluid velocity is given by 

(15) u(y) = -*E. 

ft 

It is a constant profile, independent of position along y-direction in the 
nozzle, as the yellow arrows shown in Fig. [7^b). 

It would be much easier to calculate the formation time as velocity is 
constant, since no integration is needed. Now, substituting Eqs. (5) and 
(15) into Eq. (8), we can get the formation time of a droplet as 

nR s V y droplet _ 27/xd 2 Singe 

[ j f ° rm ~ uA ~ zA E «a? - aU 2 Pe et " 

Eq. (16) suggests the formation time of a droplet is a function of surface 
tension, viscosity, gap between electrodes, contact angle, nozzle size, applied 
voltage, permittivity, etc. 



6 Simulations and Results 



Here, nozzle size a = 10//m, applied voltage U = 1500V, gap d = lOO^m. 
So, E = 1.5 x 10 7 V/m. The liquid used in the simulation is water, the 
common solution to most drug solutions. Its surface tension 7 = 72 x 
10~ 3 N/m, and after Teflon coating, the contact angle increases, 6 C = 110°. 
Assume the spatial charge density of a droplet . Simulations of an electric 
field distribution and a droplet formation process in this EHD based system 
are shown in Fig. [8} 



10 



EHD BASED PULSATILE DRUG DELIVERY SYSTEM 



(a) 



(b) 




Figure 8. 

COMSOL simulation of electric field distribution between metal needle 
and ground electrode, (a) numerical analysis domain, (b) initial setup, (c) 
beginning of droplet formation process, (d) a droplet formed by EHD. 



We also simulated the neutralization process, which is necessary for a 
newly formed droplet before ejecting into the body. Highly charged droplets 
formed by EHD method and the neutralization region with opposite charges 
are shown in Fig. [9} 

Different ratios of the spatial charge density on droplets and neutralization 
region, affect the neutralization effect, as shown in Fig. [TUJ The ratio ranges 
from 4000:1 to 1:2. Neutralization works well when the ratio is around 200:1, 
at which a formed droplet is neutralized completely. Otherwise, the charges 
generated by EHD on droplet are either incompletely neutralized or changed 
to the opposite because of the higher charge density in neutralization region. 
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(a) (b) 




Figure 9. 

Simulation of the neutralization region, (a) high charged droplets without 
neutralization, (b) neutralization region with opposite charges. 



4000:1 2000:1 200:1 




Figure 10. 

Effect of neutralization due to ratios of spatial charge density on droplets 

and neutralization region 
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Table 1. Dimensions and system data used in calculation 



Dimensions of device (|im) 



O.D. 


30 


Wpdms 


2 


ID. 


20 


Lpdms 


3 


Lueedle 


150 


a 


10 




400 


d 


100 


External voltage (V) 


1500 


Drug supply rate (kg/ min) 


0.02 


Spatial charge density (C/m 3 ) 


0.1 


Surface electric potential (V) 


0.1 



7 Discussion 

This EHD based drug delivery system can deliver a single droplet as small 
as 2.83nL in 8.5ms with a total device of 2 x 2 x 3mm 3 , and a working 
voltage of 1500^. The droplet generation is repeatable and controllable. 

The dimensions of device and data of the drug delivery system used in 
calculation are shown in Table 1. Except for DI water used in calculation and 
simulation, many other solutions could be used for this design. Properties 
of various solutions can be found in Ref421j. 

Formation time of a droplet, in Eq. (16), consists of the applied voltage, 
gap between tip of nozzle and ground electrode, nozzle size, contact angle, 
and some fluid properties, such as viscosity, permittivity and surface electric 
potential. For a perticular drug fluid, its properties are fixed. Therefore, 
the factors to change the formation time are dimensions of the system. 

Figure 11 shows the different effects on formation time caused by the four 
easily manipulated parameters in experiments, external supplied voltage U, 
gap between two electrodes d, nozzle size a and contact angle 9 C . In Fig. 
11(a) and 11(c), tf orm reduces as U and a increases separately. It means that 
it is easier and faster to form a droplet due to EHD at a higher voltage or a 
smaller nozzle size. In Fig. 11(b), tf orm increases when d grows, resulting in 
a decreasing electric field strength E. Contact angle increases after Teflon 
coating. In Fig. 11(d), only when 6 C > 90°, it can reduce tf orm . 

8 Conclusion 

Electrohydrodynamic generation, a commonly used method in BioMEMS, 
plays a significant role in the pulsed-release drug delivery system for a 
decade. 
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nozzle size, a ipm) contact angle, c (degree) 

Figure 11. 

Effects on formation time i/ orm caused by external supplied voltage U, gap 
between electrodes d, nozzle size a and contact angle 6 C 



In this paper, an EHD based drug delivery system is well designed for 
the stable, repeatable and controllable generation of a nano/micro droplet 
at a constant rate. Theoretically, we derive the expressions for the size and 
the formation time of a droplet generated by EHD method, while taking 
into account the mass supply rate and properties of drug, gap between two 
electrodes, nozzle size, and charged droplet neutralization. 

Even with an inserted metal tube or needle to enhance the efficiency 
and stability of the droplet formation, its fabrication and critical alignment 
remain challenging. 
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